Is LaOi-^F^FeAs an electron-phonon superconductor ? 
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In this paper we calculate the electron-phonon coupling of the newly-discovered superconduc- 
tor LaOi-^F^FeAs from first-principles, using Density Functional Perturbation Theory. For pure 
LaOFeAs, the calculated electron-phonon coupling constant A = 0.21 and logarithmic-averaged fre- 
quency u>i n — 206K, give a maximum T c of 0.8 K, using the standard Migdal-Eliashberg theory. For 
the F— doped compounds, we predict even smaller coupling constants, due to the strong suppression 
of the electronic Density of States at the Fermi level. To reproduce the experimental T c = 26K, 
a 5-6 times larger coupling constant would be needed. Our results indicate that electron-phonon 
coupling is not sufficient to explain superconductivity in the newly-discovered LaOi-zF^FeAs su- 
perconductor, probably due to the importance of strong correlation effects. 

PACS numbers: 71.38.-k, 74.25. Jb, 74.25.Kc, 74.70. Dd 
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The very recent report of superconductivity with the 
remarkable T c of 2&K in LafOi-^F^jFeAs [l] has stim- 
ulated an intense experimental and theoretical activity, 
aimed at identifying the possible superconducting mecha- 
nism. This compound belongs to a family of quaternary 
oxypnictides of the form LnOMPn, where Ln=La and 
Pr, M=Mn, Fe, Co and Ni; Pn=P and As, synthesized in 
1995. [2| Superconductivity was first reported in LaOFeP, 
with a relatively low T c of ~ 7 K Q , and later in F-doped 
LaOFeAs, with a maximum T c of 26 K at x=0.12 (ap- 
parently pure LaOFeAs shows no superconductivity). 

The first bulk measurements on a sample with x = 0.1 
have shown that F-doped LaOFeAs has a relatively small 
in-plane coherence length (£ a b = 81 A) and a T-dependent 
Hall coefficient [3], the electronic specific heat displays 
a vanishingly small jump at T c , and its behavior under 
magnetic field |f| as well as point-contact spectroscopy 
suggest the presence of nodes in the superconducting 
gap. All these observation suggest a strong analogy with 
the high-T c cuprates. 

A recent LSDA calculation predicts that pure 
LaOFeAs is on the verge of a ferromagnetic instability, 
due to a very high Density of States (DOS) of Fe d elec- 
trons at the Fermi level. [H A DMFT calculation in the 
paramagnetic phase, including strong-correlation effects 
beyond LDA, shows that for [7 = 4 eV, a large amount 
of spectral weight is shifted away from the Fermi level, 
and the undoped system has a bad metallic behaviour. Q 
Both papers rule out standard e — ph theory as a possi- 
ble explanation for superconductivity, without estimat- 
ing the magnitude of the e — ph coupling constant. 

In this Letter, we calculate from first-principles the 
electron-phonon properties of LaOFeAs, using Density 
Functional Perturbation Theory [13, Similar cal- 

culations, in conjunction with Migdal-Eliashberg theory, 
reproduced the superconducting properties of many stan- 
dard e — ph superconductors llj, including MgB2 12 1 



with considerable accuracy. On the other hand, they fail 
dramatically in the the High-T c cuprates, ll4H where 
the Local Density Approximation is not sufficient to de- 
scribe the strong local electronic correlations, and their 
interaction with phonons. [Pol ] 

Our calculations show that LaOFeAs is intrinsically 
a very poor e — ph superconductor, with a very weak 
e — ph coupling distributed evenly over several phonon 
branches. For electron-doped LaOFeAs we calculate an 
upper limit for the e — ph coupling constant A ~ 0.21, 
which, together with uji n — 206X, is a factor 5 too small 
to account for the observed T c = 2QK . 




FIG. 1: (color online) Crystal structure of LaOFeAs. 

LaOFeAs crystallizes in a tetragonal crystal struc- 
ture (space group 129), with a = 4.035(3.996)1, 
c = 8.741 (8. 636) A; [l6| La and As atoms occupy 2c 
Wyckoff positions, with z = 0.1415(0.1413) and z = 
0.6512(0.6415) respectively; O and Fe atoms occupy 2a 
and 2b Wyckoff positions. 

The structure, depicted in Fig. [TJ consists of alternat- 
ing Fe-As and La-0 layers. Fe and O atoms sit at the 
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center of slightly distorted As and La tetrahedra; the 
As tetrahedra are squeezed in the z direction; the Fe-As 
distance is 2.41 (2.34) A, and the As-Fe-As angles are 
either 107.5 (105.81) or 113.5 (117.1) degrees. Fe atoms 
also bond to other Fe atoms in plane, which are arranged 
on a square lattice at a distance of 2.85(2.83)A. 

The gross features of the band structure of LaOFeAs 
are very similar to those of LaOFeP. 0, B 0] • Measur- 
ing energies from the Fermi level, O p and As p states 
form a group of 12 bands extending from ~ —6 to —2 
eV. La-/ states are found at higher energies, at ~ 2 eV. 
Apart from a weak hybridization of the ti g states with 
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FIG. 2: (color online) Band structure of LaOFeAs, decorated 
with partial characters of the e g (top) and tig (bottom) Fe-d 
bands. The orientation of the coordinate system is chosen 
so that Fe-Fe bonds are directed along x, y axes; the zero of 
the energy coincides with the Fermi level. The arrows indicate 
the splitting induced by the elongation/shrinking of the Fe-As 
tetrahedra (see text). 

As p states, at —3 eV, the ten Fe— d states are localized 
in an energy window extending ±2 eV around the Fermi 
level, where they give the dominant contribution to the 
DOS. The derived bands do not split simply into a lower 
e g (d x 2_ y 2 and d 3z 2_ 1 ) and higher ti g manifold, as pre- 
dicted by crystal field theory (see Fig. Due to the 
presence of Fe-Fe directed bonds, d x 2_ y 2 orbitals, which 
lie along the bonds, create a pair of bonding-antibonding 
bands located at —2 and +1 eV. d 3z 2_ 1 bands states split 
into two sub-bands. Due to the distortion of the Fe-As 
tetrahedra, d xy states are inequivalent to d xz , d yz states, 
and the resulting the ti g bands form a complicated struc- 
ture centered at ~ —0.5 eV. 

The Fermi level cuts the band structure in a region 
where the DOS is high ( 2.1 states/eV spin) and rapidly 
decreasing; a pseudo-gap opens in the electronic spec- 
trum around 0.2 eV. As pointed out in previous publi- 



cations such a high DOS at the Fermi level drives the 
system close to a magnetic instability. 0, Q 

The Fermi surface comprises a doubly-degenerate 
cylindrical hole pocket centered at the T point, and a 
doubly-degenerate electron pocket centered at the M 
point; these sheets have a dominant d xz ,d yz character. 
A small 3D pocket centered around the F point is also 
present (see Fig. 3 of Ref. 3). The plasma frequencies 
are strongly anisotropic (uj xx = 2.30, uj zz = 0.32 eV). The 
distortion (elongation or shrinking) of the Fe-As tetrahe- 
dra modulates the splitting of the two d 3z 2_ 1 bands, and 
the relative splitting between xy and xz, yz bands along 
the T-Z line, as indicated by the small arrows Fig. [2] 
A 1% percent compression of the tetrahedra along the 
c direction changes the splitting of the two d 3z 2_ 1 by 
~ O.ley, driving one of them closer or further from the 
Fermi level. This explains why the position of the idz 2 — 1 
band, and the weight of the associated DOS, varies in lit- 
erature, depending on the crystal structure used. [1, Q 

Fig. [3] summarizes the e — ph properties of LaOFeAs; 
the results refer to pure LaOFeAs in the paramagnetic 
phase. It has been shown that the pure compound is 
close to a magnetic instability and to a metal-insulator 
transition due to electronic correlation. 0, Q Electron 
doping strongly suppresses the tendency to magnetism 
and reduces strong correlation effects, and assuming a 
paramagnetic ground state is probably appropriate for 
the F-doped compound. Also, we checked by calculations 
in the virtual crystal approximation that the effect of F 
doping is well described by a rigid-band model, and the 
only effect of doping is a rigid-band shift of the Fermi 
level, in a region where the electronic DOS is lower (a 
10% doping correspons to a 40 % reduction of the DOS). 
Therefore, the results for the undoped compound can 
be considered representative also for the electron-doped 
compound, provided that the reduced DOS is taken into 
account. 

In the left panel of Fig. [3l we show the calculated 
phonon dispersion relations of LaOFeAs; our frequencies 
are in very good agreement with those of Ref. § , where a 
slightly different crystal structure was used. In the mid- 
dle panel of the same figure, we show the atom-projected 
phonon DOS. The spectrum extends up to 500 cm -1 ; the 
vibrations of O atoms are well separated in energy from 
those of other atomic species, lying at ui > 300cm -1 . 
The vibrations of La, Fe and As occupy the same en- 
ergy range, and the eigenvectors have a strongly mixed 
character. Similarly to the electronic bands, the phonon 
branches have very little dispersion in the z direction. 
Analyzing the evolution of the phonon eigenvectors in 
the Brillouin Zone (BZ) reveals that there is no clear 
separation between in and out-of-plane vibrations, as it 
often happens in layered compounds. The three major 
peaks in the phonon DOS at u> = 100, 200 and 300 cm -1 
do not show a definite in-plane or out-of-plane character, 
and cannot be easily traced back to a single vibration pat- 
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tern. This complicates the interpretation of the Eliash- 
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FIG. 3: (color online) Electron-phonon properties of 
LaOFcAs. Left: Phonon dispersion relations; the radius of 
the symbols is proportional to the partial A of each phonon 
mode, Eq. ((3J. For T and Z points, where the coupling di- 
verges numerically, we use a different scaling factor. Middle: 
Atom-projected phonon DOS. The projection on in- and out- 
of-plane modes (not shown) does not show any clear sepa- 
ration between the patterns of vibration. Right: Eliashberg 
function a 2 F(uj) (full line) and frequency- dependent electron- 
phonon coupling A(w) (dashed line), Eqs. [TJ2] 

berg spectral function a 2 F(u>), shown in the rightmost 
panel of Fig. [31 together with the frequency-dependent 
e—ph coupling constant A(u>): 
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A comparison of the Eliashberg function with the phonon 
DOS shows that, except for the high-lying O modes, 
which show very little coupling to electrons, the e — ph 
coupling is evenly distributed among all the phonon 
branches. Low- frequency phonons around 100 cm -1 pro- 
vide ~ 75% of the total A, due to the 1/fl factor in Eq.[2j 
but the e — ph matrix elements g are comparable for all 
group of phonons. 

It is interesting to note that this almost perfect propor- 
tionality between the Eliashberg function and the phonon 
DOS is never encountered in good e—ph superconductors, 
where the coupling to electrons is usually concentrated 
in a few selected phonon modes. This is best explained 
in terms of phonon patterns that awake dormant e — ph 
interaction between strongly directed orbitals. 

An extreme example in this sense is MgB2, which 
achieves a T c of 39 K thanks to a strong coupling be- 
tween bond-stretching phonons and strongly covalent a 



bands, but the same applies also to more traditional su- 
perconductors, such as the A15, NbC, and even normal 
metals. 

In LaOFeAs, all phonon modes give a comparable, 
small contribution to the total A; this indicates that there 
are no patterns of vibration with a dramatic effect on the 
electronic band structure around the Fermi level. A pos- 
teriori, this is not surprising since in LaOFeAs the only 
bands derived from directed bonds, (d x 2_ y 2 in Fig. [T}, 
which could experience strong coupling to Fe vibrations 
in plane, sit far from the Fermi level. 

The distribution of the coupling is also shown in the 
left panel of Fig. [3j where the radius of the circles is pro- 
portional to the mode A, i.e. to the partial contribution 
of each phonon mode to the total e-ph coupling 
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where 7„ q are the e-ph linewidths; summing Eq. |[3} over 
the phonon branches v and averaging on the BZ give the 
total e-ph coupling A. The circles are evenly distributed 
over several phonon branches. The largest couplings are 
concentrated around the T(Z) points, where the intra- 
band nesting is large, and around the M point, where 
the interband nestig between the hole and electron cylin- 
ders take place. 

The total e—ph coupling constant A, obtained by nu- 
merical integration of Eq. ([2]) up to uj = oo, is 0.21; 
this, together with a logarithmically-averaged frquency 
ujin = 205K, and fi* = 0, gives T=0.5 K as an upper 
bound for T c , using the Allen-Dynes formula fl7j |. Nu- 
merical solution of the Eliashberg equations with the cal- 
culated a 2 F(ui) function gives T c =0.8 K. To reproduce 
the experimental T c = 26K, a five times larger A would 
be needed, even for fi* = 0. Such a large disagreement 
clearly indicates that standard e—ph theory cannot be 
applied in LaOFeAs, in line with recent theoretical works 
which emphasize the role of strong electronic correlations 
and/or spin fluctuations. [1, Q 

The numerical uncertainity on the calculated value of 
A, connected to limited sampling of the BZ in k (elec- 
trons) or q (phonons) space integration, is at most 0.1, 
and definitely not sufficient to raise A to ~ 1.0. We fur- 
ther notice that electron doping, reducing the DOS at 
the Fermi level, without introducing new bands at E^, 
would further reduce the value of A. Therefore, the value 
A = 0.21 for the undoped material is actually an upper 
bound for the value in the e-doped compound. This value 
is lower than what is encountered in any known e — ph 
superconductor; for comparison, A = 0.44 in metallic alu- 
minum, where T c = 1.3 K. 

In LaOFeAs, both the electronic DOS at the Fermi 
level, and the value of the average phonon force constant 
are in line with those of other e—ph superconductors. 
The occurrence of a small A is due to its extremely small 
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matrix e — ph elements, connected to the strongly delo- 
calized character of the Fe-d states at ±2 eV around the 
Fermi level. This is an intrinsic property of this material, 
which could hardly be modified by external parameters, 
such as pressure or doping. For the same reason, our re- 
sult is quite stable with respect to the minor differences 
in the electronic structure around the Fermi level, which 
have been observed in literature. 

In principle, multiband and/or anisotropic coupling 
could provide the missing factor 5 in the coupling, but 
this is very unlikely to occur because it would require a 
very large anisotropy of the distribution fl8i | . Other in- 
teractions, repulsive in the s-wave channel but attractive 
in the d- or p-wave one, may increase T c . 

In conclusion, we have calculated the electron- 
phonon properties of the newly discovered superconduc- 
tor La[Oi_ a; F a ;]FeAs using Density Functional Perturba- 
tion Theory. The undoped compound is close to a mag- 
netic instability, due to the presence of a very sharp peak 
in the electronic Density of States. Doping with electrons 
moves the system away from the magnetic instability, re- 
ducing the DOS at the Fermi level, without altering the 
band structure substantially . 

Despite the high value of the DOS at the Fermi level, 
the calculated value of the e-ph coupling constant for the 
pure compound is only A = 0.21, which is a factor 5 too 
small to yield the experimentally measured T c within the 
scope of standard ME theory. 

The value A = 0.21 in LaOFeAs is very close to those 
which have been estimated using LDA in the supercon- 
ducting cuprates (A ~ 0.3). Q3] Similarly to the 
superconducting cuprates, in this compound supercon- 
ductivity cannot be described using standard LDA cal- 
culations and Migdal-Eliashberg theory. Although our 
findings do not necessarily imply that superconductivity 
in e-doped LaOFeAs is due to an exotic mechanism, they 
clearly indicate that strong correlation effects beyond the 
LDA play an important role and must be included in any 
realistic description of this material. 

Technical Details: For the atom-projected band and 
DOS plots in Fig.[2]we employed the full-potential LAPW 
method [l9[ as implemented in the Wien2k code. [2(| 
Calculations of phonon spectra, electron-phonon cou- 
pling and structural relaxations were performed using 
planewaves and pseudopotentials with QUANTUM- 
l|, [22j. Whenever possible, we cross-checked 
given by the two codes and found them 
to be in close agreement; for consistency, we used the 
same GGA-PBE exchange-correlation potential in both 
cases. [IJ] 
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